Current and noise in a model of an AC-STM molecule-metal junction 
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The transport properties of a simple model for a finite level structure (a molecule or a dot) 
connected to metal electrodes in an alternating current scanning tunneling microscope (AC-STM) 
configuration is studied. The finite level structure is assumed to have strong binding properties 
with the metallic substrate, and the bias between the STM tip and the hybrid metal-molecule 
interface has both an AC and a DC component. The finite frequency current response and the 
zero frequency photo-assisted shot noise are computed using the Keldysh technique, and examples 
for a single site molecule (a quantum dot) and for a two-site molecule are examined. The model 
may be useful for the interpretation of recent experiments using an AC-STM for the study of both 
conducting and insulating surfaces, where the third harmonic component of the current is measured. 
The zero frequency photo-assisted shot noise serves as a useful diagnosis for analyzing the energy 
level structure of the molecule. The present work motivates the need for further analysis of current 
fluctuations in electronic molecular transport. 



I. INTRODUCTION 



The understanding of electronic transport phenom- 
ena in molecule-electrode interfaces is a key element of 
any technological realization of the potential of single- 
molecule electronics [1,2]. The basic physics governing 
the current voltage characteristics in the coherent regime 
is reasonably well understood in terms of generalizations 
of the Landauer model where current is viewed as a quan- 
tum scattering problem and the effect of the electrodes 
is included via a self-energy contribution. Such a model 
has been successfully used to interpret data from both 
isolated molecules in STM, AFM and break junctions ex- 
periments [3-5] , as well as experiments on Self Assembled 
Monolayers [6-8]. Most experiments on molecular-scale 
devices have been carried out on a two-terminal assem- 
bly, but measurements have also been reported in gated 
systems [9,10]. A vast majority of experimental and the- 
oretical studies have focused on the measurement of the 
current in such molecular junctions. Generally speak- 
ing, in the coherent tunneling regime, the zero frequency 
current 7(0) in molecular wire junctions depends on the 
electronic structure of the extended moZecufe-electrode 
system, the spatial profile of the voltage drop, the electro- 
chemical potential across the junction and the self-energy 
of the contacts that contains the information about the 
coupling between the finite molecule and the extended 
electrode [11,12]. Yet from the point of view of meso- 
scopic physics, the current alone is not sufficient to fully 
characterize the transport: the second moment of the 
current - the noise - has been shown to provide crucial 
information concerning the effective charge of the carri- 
ers and their statistics [13,14]. While noise measurements 
still represent a certain challenge in molecular transport - 
as experiments are typically performed with tunnel con- 
tacts - recently noise was computed for simple molecules 



using a Hiickel type Hamiltonian [15,16] and in the case 
of Coulomb interaction in the molecule using an Ander- 
son impurity model [17]. The present work deals with 
the computation of both current and noise when the bias 
imposed across the molecule has both a constant (DC) 
and an alternating (AC) component [18,19]. 

Indeed, the alternating current STM (AC-STM) was 
suggested by Kochanski [20] as a experimental technique 
which could extend the capabilities of conventional STM 
from conducting to insulating surfaces. A recent exper- 
iment [21] on the oxidation of lead sulfide surfaces con- 
firmed that crucial information could be obtained from 
the analysis of the third harmonic of the current. Never- 
theless, there has been so far no report of an AC-STM 
diagnosis on a Self- Assembled Monolayers as a means to 
probe single molecule transport. 

From the point of view of the theoretical community, 
the problem of electron transport in the presence of both 
a DC and an AC bias voltage is often referred to as photo- 
assisted transport [22-25]. In the absence of an AC bias, 
most treatments of molecular wire junctions focus on the 
problem of stationary current and conductance. On the 
other hand, the problem of AC-STM requires an explicit 
time-dependent approach that allows a full description of 
the transient terms, as well as the stationary terms. In 
this regard a Hamiltonian formulation of transport such 
as the Keldysh method [26-28] is more powerful than 
the scattering-based Landauer approach [29] . It allows a 
proper treatment of the electron statistics in real time. 

Here a simple model junction is considered, which con- 
sists of two metal electrodes and a finite electronic struc- 
ture, either a dot or a molecule described by a tight- 
binding Hamiltonian, which is coupled to the electrodes 
through a single contact. The system is subject to the 
combined action of a DC voltage bias and a time- varying 
field. The Keldysh formulation allows to compute both 
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the frequency-dependent current response as well as the 
photo-assisted shot noise. The molecule is assumed to 
be well connected to one of the electrodes (the substrate) 
and its binding properties are computed non perturba- 
tively. The coupling to the other electrode (STM tip) is 
treated to second order in the hopping amplitude, but 
the field amplitude appears to all orders in this problem 
as will be shown below. Under simplifying assumptions 
(Poisson regime), the noise is just proportional to the av- 
erage stationary current itself, and therefore both contain 
essentially the same information. This Shottky relation 
between current and noise ceases to be valid for higher 
order harmonics of the current and therefore the study 
of the noise can reveal important informations about a 
periodically driven system. Finally, one should empha- 
size that the general method employed here is in no sense 
specific to a mctal/single-dot/metal junction. Here one 
has chosen to proceed with specific computations on this 
particular system, but additional results on the two site 
molecule can be generalized to any linear molecular chain. 

This article is organized as follows: Section 2 presents 
a description of the model. Section 3 is devoted to the 
computation of non-equilibrium transport quantities. In 
section 4, the main results for the photo-assisted shot 
noise for both the case of a single dot molecule and the 
two-site molecule are shown and discussed. The higher 
harmonics of the current are examined in section 5. Fi- 
nally, section 6 contains some concluding remarks and 
speculations. 



II. MODEL HAMILTONIAN 




FIG. 1. System : the molecule (quantum dot) sandwiched 
between a substrate and an STM tip. 



Consider the model for a metal-moZecit/e-metal junc- 
tion depicted in Fig. 1. The substrate and the tip are 
assumed to be metallic and are both described by an infi- 
nite bandwidth model. Both are therefore characterized 
by constant density of states ps and px- 

In turn, the molecule will consist either of a single dot 
or of a conjugated structure (described by a tight bind- 
ing Hamiltonian) . It contains sites numbered from 1 to 
D. The connection between the molecule and the sub- 
strate occurs on site 1, while the connection between the 
molecule and the tip is located at site D. This general 
approach allows to compute the transport properties for 
a molecule with an arbitrary number of sites, but specific 
plots will be obtained for D = 1 and D = 2. 

To be more specific, the connection between the 
molecule and the substrate is assumed to be quite strong. 
Indeed, such strong binding to a substrate is known to 
occur when the molecule is functionalized: when for in- 
stance an end group of atoms in the molecule is replaced 
by a sulphur atom, which has strong binding properties 
on a gold substrate. 

When the other electrode (the STM tip) is brought 
close to the molecule, and a bias voltage is imposed be- 
tween both metallic electrodes, non equilibrium flow of 
electrons occurs. The key point of this work is to use an 
alternating current STM for probing transport. This is 
modeled by adding an alternating contribution to the DC 
bias Vo imposed on the junction: V(t) = Vq + V\cos(ujt) 
where V\ is the amplitude of the AC bias voltage and lu 
is its frequency. 

The sites involved in the junction between the molecule 
and the tip are denoted D (as in dot) and T (tip). With 
this convention, the Hamiltonian describing tunneling 
through this junction reads: 

H TD (t) = T(t)c f T c D + T*(t)c f D c T (1) 

where c/c^ are fermionic annihilation/creation opera- 
tors. The time dependence of the voltage bias is thus 
included in the hopping term using the Peierls substi- 
tution: starting from a formulation where the electric 
field at the junction is described solely by a scalar po- 
tential, a gauge transformation with gauge parameter 
e x(t) c = — e J V(t)dt is performed. The phase of the 
amplitude has a linear dependence associated with the 
DC bias, but it also contains a modulation due to the 
AC voltage: 

r(t) =T° exp[-ie X (t)/c] 

= r£ cxp(i (u t + ^-sm(wi)j ) (2) 

where r!^ is the time independent tunneling amplitude 
of the hopping between the tip and the molecule, and 

<*>o,i = eV o,i- 

An advantage of the Peierls substitution is that the 
current operator can be directly obtained by differentiat- 
ing the hopping Hamiltonian with respect to the gauge 
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parameter: 



I(t) = -c 



dH TD (t) 
dx(t) 



(3) 



Moreover, the current-current correlations (the cur- 
rent noise) are defined by 

S{t, t>) = \ [(J(t)J(f ) + J(f )/(*)> - </(t)>(J(f )>] (4) 

Note that this correlator is symmetrized, which is appro- 
priate when computing zero frequency noise: when con- 
sidering the finite frequency spectral density of noise, this 
ceases to be true and one needs to consider the unsym- 
metrized correlator, which makes the difference between 
emission or absorption [22]. 



III. NON-EQUILIBRIUM TRANSPORT 

Note that because of the time dependent bias, any 
straightforward extension of the Landauer approach has 
to be ruled out, and the Kcldysh formalism is employed. 

The expectation value of the current (Eq. (3)) can 
be directly cast in terms of the spectral Green's func- 
tions. For the noise (Eq. (4)), we can factorize the two 
body correlation function following a mean field decou- 
pling scheme, which is exact for non-interacting electrons. 
Green's functions are represented by 2 x 2 matrices: each 
time argument can be pinned on either branch of the 
Kcldysh time contour. However, recall that these Green's 
functions are not independent. For convenience, here, 
one chooses to use the version of Dyson equations for 
advanced, retarded and spectral Green's functions : 

G a a0 (t,t/)=iQ(t?-t)({c a (t),c+(t>)}) (5) 

G r af3 (t,t/) = -ie(t-t')({c a (t),c+(t')}) (6) 

G+p(t,t')=i(cl(t>)c a (t)) (7) 

G-+(t,t') = -i(c a (t)cl(t>)) (8) 

Here a, (3 are terminal indices. Inserting these Green's 
functions in the current and the real time current-current 
correlator yields: 

J(t) = 2eRe [T(t)G+ T (t, *)] 
S(t,t') = 2e 2 Re\r(t)T*(t r )G+ T (t',t)G D +(t,t')) 

+ T*(t)T(t')G+ D (t',t)G T +(t,t')] (9) 

Here one is only interested in the lowest order contribu- 
tion in the tip-molecule hopping amplitude. The Dyson 
equation for the above Green's function reads (omitting 
terminal indices summations and time integrations): 

G a,r = g a,r^ + ^TQa.r^ 

G + -'- + = (1 + G r Y, r )g + -'- + {l + £ a G Q ) (10) 



In the perturbative scheme, the G's are the dressed 
Green's functions -which have both diagonal and non- 
diagonal non-zero elements in the terminal indices- while 
the g's are the bare Green's functions with only di- 
agonal non-vanishing elements. The self-energies £ 
are related to the hopping amplitude : S r = T(t), 
S a = T*(t), S+--+ = 0. At this point the time- 
dependent hopping amplitude is expanded using the gen- 
erating function of the Bessel functions [30] T(t) = 

r Y^+oc i(u +nw)t j ( un\ 

1 T Z^n=-oo e J ™ V w )■ 

The strong coupling between the substrate and the 
molecule allows to define a compound moZecwZ^substrate 
system. The Green's function of the molecule becomes 
dressed by the substrate-moZecwZe coupling, which can 
be treated non perturbatively. In the following, G de- 
notes the non-equilibrium dressed Green's functions for 
this composite system. 

A Fourier transform (h(uj) = dte tbJt h(t)) helps to 
characterize the current response at finite frequency. Be- 
cause time-translational invariance is broken by the AC 
bias, a double Fourier transform is implied for the noise: 



n,m ~ 00 

xj [5(fl + (n-m)uj) + 8(fl-(n-m)u;)} 

x G^(w 2 )5tt( w 2 -uiQ-mw) 
+ 6 (CI + {n- m)uj) 

x G r DD (Lo 2 )gTT ( w 2 - w - raw) 

- G a DD (ijj-2)gTT ( w 2 - w - nw)] I 

/+oc 
dcj 2 
-co 

g^xi^ + 0' — w — rauS) 

+ g T ~T~(u>2 — 0' — djQ — nuj) 
g^xi^z + ^ — w o — niuj) 

+ g T : T ~(oj2 — O, — lu — nuj) j (11) 



At this point, it is convenient to introduce the tip's 
Green's functions (Appendix A 1) where pr is the den- 
sity of states in the tip. The noise and the current thus 
become: 
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m = le7 rp T \r T \^Jn(^) J m (^) 
x I [5(0, + (n- m)w) + 5(0 - (n - m)w)] 

/+oo 
du 2 G DD (u) 2 ) 
-oo 



+2 5(fi + (n-m)w) 



^2 (G5, D (W2) - GS, D (W2) 



/ (Lj 2 G DD (u 2 ) \ (12) 



S(Cl,Cl') = (2 e ) le 7rp T |r^| 2 ^ J„ (^) J m (^) 

n.m 

x5(Cl + Cl' + (n- m)u) 

r /" + °° 

x j -2 duj 2 G dd (uj 2 ) 

+2 / du 2 G D +(u 2 ) 



dcj 2 



ujQ+mu) — H 



(13) 



The noise and current are fully determined once the 
dressed Green's functions of the molecule (dot) are spec- 
ified. Note that Eqs. (12) and (13) are valid for an 
arbitrary model of the mo/ecw/e/substrate, which could 
in principle include interactions. The S functions in Eq. 
(12) and (13) restrict the frequency dependance of the 
current to integer multiple of the driving frequency u>. 
The mth harmonic of the current, I(mu>), is obtained 
by integrating Eq. (12) on the detector resolution s : 
I(muj) — f™"*g/ 2 dill (CI), then getting rid of the S func- 
tion singularities. The detector resolution s is of course 
much smaller than all the frequencies of this problem. 

Here we focus on a single dot and on a two site 
molecule. The dressed Green's functions are derived in 
Appendices B 1 and B 2. For the single dot they read: 



G a 6n(Cl) = 



1 



G+- D (Cl) = 2m Ps T% 



f(Cl - E F ) 



(CI - E D ) 2 + [ixpsYlf 



G DD (Cl) = -2mp s T 2 s 



(l-f(Cl-E F )) 

(ci - E D y + (7T Ps r|)2 



(14) 



The coupling to the metallic substrate imposes a width 
TrTgps- Furthermore, the occupation of the dot (as it ap- 
pears in G DD ) is specified by the Fermi function / of the 
substrate. 



For two site molecule, the dressed Green's functions 
are derived in Appendix B 2 : 

G a J D (Cl) = - n -^ T s 



G DD (Cl) = 2z7rp s r 



t 2 - Cl(Cl - iTrp s T 2 s ) 

t 2 f(Cl - E F ) 



s \t 2 -Cl(Cl + iirp s r 2 )\ 2 



G DD (Cl) = -2znp s T 2 s 



2 t 2 [l-f(Cl-E F )} 



\t 2 -Cl(Cl + iirp s r 2 )\ 2 



(15) 



where t is the hopping between the two sites of the 
molecule. 

Note that, in both cases (single dot/ double dot), the 
spectral dressed Green functions may be written in terms 
of the density of states of the dressed system : 



G DD (Cl) = 2iIm[G a DD (Cl)]f(Cl - E F ) 
G D +(Cl) = -2iIm[G a DD (Cl)][l - f(Cl - E F ) 



(16) 



where Im[G DD (Cl)]/'K is the density of states on site 
D . This density should be computed directly from the 
expression of the advanced dressed Green's function : 



Im[G a DD (Cl)}\ ldot 



Im[G DD (Cl)]\ 2dots = 



(CI - E D ) 2 + (np s r 2 ) 2 

irp s r 2 s t 2 
(t 2 - CI 2 ) 2 + (np s T 2 ) 2 Cl 2 



(17) 



In the following, the integrated version of the density 
of states denoted by F is needed : 



F(a,b)= j du Im[G a DD (u)] . (18) 

J a 



IV. PHOTO-ASSISTED SHOT NOISE 

One is now in a position to compute the current har- 
monics and the photo-assisted shot noise using the ex- 
plicit expressions for the dressed Green's functions. Be- 
fore displaying these results, as a reference it is useful 
to recall what happens in a tunnel junction between two 
metallic electrodes (no molecule). This allows to relate 
to the results of Ref. [31], where photo- assisted shot noise 
was studied in a mesoscopic conductor from the point of 
view of scattering theory. For the current and the noise 
our perturbation scheme yields at zero frequencies : 



I(Cl = 0) = 4e7r 2 p T p s T 2 £ J 2 (^) (uj + nuj - E F ) 

n 

5(0, 0) = (2e)Ae^p TPs T 2 s ]T J 2 (^) |wo +nuj-E F \ 

n 

(19) 
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The noise is continuous, but its derivative has jumps at 
integer values of w /w, precisely the result of Ref. [31] 
when one assumes that the transmission coefficient used 
for scattering theory does not depend on energy. It is 
therefore relevant to generalize these results to a situation 
where the "mesoscopic system" has an internal discrete 
level structure, which is the case for molecular electronic 
transport. This is expected to have repercussions on both 
the current harmonics and the noise. 

Using the general formulation of the spectral Green's 
functions (Eq. 16), the zero frequencies noise read : 

S(0,0) = 8e 2 7rp T |r° 4 (^) s W^o + nw - E F ) 



x F(E F ,w + nw) 



(20) 



where F(a, b) is defined in Eq. (18). 

In order to make a diagnosis on the spectrum, it is use- 
ful to display also the first derivative of the noise with 
respect to the DC voltage: 

= 8e 2 7rp T |rO | 2 £ \jl (^) S gn(w Q + nw E F ) 



x Im[G a DD {w + nu>)} 



(21) 



o 

to 
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— r s = o.i 

r s = i.o 

__. r s = io. 




Ef "' w /w 

FIG. 2. Zero-frequencies noise S(0, 0) (upper graph) and 
its first derivative dS(0, 0)/dwo (lower graph) as a function of 
the ratio ujo/lo in the regime of a small AC bias (lu\/lu = 0.1) 
for the single dot system. Fermi energy of the substrate Ef is 
fixed to —1.4a; and energy level of the dot Ed is the origin of 
the energy. Comparisons are made between the two electrodes 
case and different values of hopping Ys in the dot-metal junc- 
tion. 



A. Single dot problem 



The quantities which are plotted on Figs 2 to 6 are 
normalized according to the prefactors which appear on 
Eqs. (22) and (23). 



One of the motivations for studying photo- assisted shot 
noise is to probe the finite frequency spectral density of 
noise using the external frequency w of the AC bias. We 
will therefore focus on the photo-assisted shot noise with 
both frequencies in Eq. (13) set to zero, because low fre- 
quency measurements (still above, say lOOKHz, to avoid 
1/f noise) are more accessible experimentally than high 
frequency ones. Using the definition of the density of 
states for a single dot system (Eq. (17)) and after per- 
forming the integration over frequencies in Eqs. (20) and 
(21), the noise at zero frequencies and its derivative read: 



S(0,0) = 8 e 2 7rp T |rO| 2 ]Tj 2 



2 (^1 



\ W 



arctan 



( {E D -E F ) \ 
\ -PsT% ) 



arctan 



sgn(ujQ + nw — E F ) 
(E D - up- nuj) 



(22) 



— - — — = 8e 7T p T p s L s \\- T \ 2_^J n \ — 



\ ui 

n 

sgn(u>o + nw — E F ) 



(wo 



nw 



E D Y + (npsTir 



(23) 



Consider first the case where the dot level lies above the 
Fermi level of the substrate. Fig. 2 displays the noise and 
its first derivative with respect to the DC bias for a small 
value of the AC bias amplitude. From the point of view 
of our analytical results, this means that very few terms 
contribute to the sum over integers in Eq. (22): n = 
dominates. In the case of large broadening, we recover 
similar results to the two electrodes system. The dot is 
metallized to the extent that transport is fully specified 
by the two Fermi distributions of the electrodes. As the 
broadening decreases, the noise acquires a step- like struc- 
ture close to the bare dot level, which gets sharper for 
lower Ts- This is associated to the fact that the current 
increases drastically at this particular DC bias. Far from 
this DC bias value, the noise saturates because the AC 
bias does not allow for the exchange of "photon" quanta 
w. The peak appearing in the derivative of the noise for 
the case of a two electrodes system is a consequence of 
the abrupt change of sign of the current. 
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Ef 1 w °/ w 
FIG. 3. Same as Fig.2 for an AC bias: wi/w = 1.0. 

When one further increases the AC bias (Fig. 3), pho- 
tonic exchanges become relevant: this is illustrated by 
the appearance of new shoulders in the noise (peaks in 
the noise derivative) for integer values of the ratio Uo/u). 
The dominant peak in the derivative is located at the dot 
level ojo = 0. Side peaks at ujq = ±u are clearly visible, 
while those at ujq = ±2w can only be identified with the 
help of the noise derivative. These peaks are smoothed 
when the substrate coupling is increased, eventually lead- 
ing to metallization as before. 





Ef luq/uj 

FIG. 4. Same as Fig.2 for an AC bias: wi/w = 10.0. 



FIG. 5. Same as Fig. 3 in the case of a resonance between 
the energy level of the dot and the Fermi energy of the sub- 
strate: E F = E D = 0. 

The behavior of this system turns out to be strongly 
dependent on the location of the dot level. In Fig. 5, 
the dot level matches the Fermi level of the substrate. 
The AC bias amplitude has been adjusted as in Fig. 3 
where essentially wo = 0, ±ui give salliant features. First, 
notice that the noise does not display steps any more. 
Qualitatively, it is understood that current flows when 
the absolute value of the bias is large compared to the the 
"photon energy" u>: in this situation there is no balance 
between carriers emmitted from the tip to the substrate, 
and vice versa. Next, the current - and thus the noise - 
must vanish when the tip and the substrate have match- 
ing Fermi energies. More interesting is the reduction of 
the noise observed when ui n = ±u. Current can flow with 
the help of virtual transitions through the dot, but for 
this AC bias where essentially only single photon tran- 
sitions contribute, there is an additional channel which 
brings about noise reduction. This explains the two side 
dips in the noise at ojq = ±a-> in Fig. 5. It is interesting 
to note that the behavior of the noise and of its deriva- 
tive is drastically modified only when the coupling to the 
substrate is small, whereas for a metallized dot the typ- 
ical behavior associated with two metallic electrodes is 
recovered. 



B. Two— site molecule 



Turning now to the case of a large AC bias (Fig. 4), 
one notices that the step like structure is preserved but 
the width of the steps is apparently doubled compared 
to Fig. 3. This is only apparent, as the noise derivative 
displays secondary peaks associated with odd photon ab- 
sorption/emission. The relative importance of odd and 
even photon peaks can be traced to a high sensitivity of 
the noise with respect to the AC bias amplitude u>i, due 
to the oscillations of the Bessel functions (see Eq. (23)). 



The expressions for the dressed Green's functions of 
the two-site molecule are inserted in the expressions for 
the zero-frequencies noise. The two-site molecule has 
a hopping amplitude t between the two sites, so the 
bonding and anti-bonding orbitals are separated by 2t. 
Note that when the Fermi level of both the tip and the 
substrate lie between these orbitals, these two level can 
mimic the HOMO (Highest Occupied Molecular Orbital) 
and the LUMO (Lowest Unoccupied Molecular Orbital) 



G 



of a more complex molecule. Here one is interested in 
understanding how the presence of the two sites modifies 
the noise with respect to the single dot case. Expressions 
for the noise and its derivative are given in Eqs. (20) 
and (21) with the corresponding density of states (see 
Eq. (17)). 



o 

cT 




o 
o 

CO 



FIG. 6. Zero-frequencies noise 5(0, 0) for a two-site 
molecule as a function of the ratio loq/lo for AC bias 
oji/cj = 5.0. Full line , F s = 0.001; dashed line, F s = 1.0; 
dashed-dotted, V s = 3.0; dotted line, T s = 100. The hop- 
ping between the two sites is t — 3cu. Vertical dashed line 
correspond to the bonding and antibonding orbitals. Plots 
are given for: a) Ef = —1.4a;; b) Ef = t (resonance with 
highest level) 

Consider first the case where the Fermi level of the 
substrate is located in between the two molecular levels 
(top part of Fig. 6). For sharp levels (small coupling 
to the substrate), the noise has a staircase structure, as 
expected from photo-assisted transport, but contrary to 
the single dot case, the steps are not monotonous. The 
noise has minima close to the location of the bonding 
and antibonding orbitals. When the levels are broad- 
ened by a better coupling to the substrate, these steps 
are smoothed out, but the minima survive at first. Fur- 
ther increasing the coupling to the substrates gives rises 
to monotonous steps which resemble the noise charac- 
teristic of a single dot with a sharp level. This can be 
explained as follows. As long as the substrate coupling 
Ts is smaller than the hopping between the two sites 
of the molecule, increasing Ts smoothes out the photo- 
assisted shot noise. But for Ts > t, the system behaves 
as if it was a single dot (the dot closest to the tip) . This 
interpretation is confirmed by comparing the density of 
states of the double dots/single dot case (Eq. (17)) : the 
limit of a very large broadening Ts >> t for the double 
dots density of states is simply the single dot density of 



states with a width t 2 /(irpsT s )(<< 1). Therefore, in- 
creasing Ts brings sharper features in the noise rather 
than smoothing them. Then, the noise associated with 
the large coupling in Fig. 6, is quite similar to the stair- 
case encountered in Fig. 4. 

Turning now to the situation where the Fermi level 
of the substrate coincides with the highest orbital of the 
molecule, one notices a noise characteristic which is again 
similar to a single dot level, except for the fact that it is 
inverted. For weak coupling to the substrate, one gets 
sharp steps, but as these steps are smoothed upon in- 
creasing the coupling, one recovers the same crossover 
back to the sharp steps of a single level. It does seem nat- 
ural that if it is now the lowest level of the molecule which 
is pinned to the substrate Fermi energy, the noise char- 
acteristic is inverted with respect to Fig. 6b (not shown) . 
In particular, for a very strong substrate coupling, we 
recover curves similar to Fig. 3. 



V. HIGHER HARMONICS OF THE CURRENT 

Using the definition of the density of states in Eq. (17), 
Eq. (12) gives for the mth harmonics of the current : 

= muj) = 2eTr PT \T° T \ 2 J n (— ) 

n 

n+m (~~) (F(-oo,Lo + (n+m)Lo)+F(E F ,LOo+nuj)) 



X < J, 



(24) 

The stationary current is given by taking m = in this 
expression : 

I(fl = 0) = 4e7T PT \T° T \ 2 Jn (— ) F{E F ,w + mo) 

n 

(25) 

According to Eq. (20), Shottky relation occurs between 
the zero frequencies noise and the stationary current : 
5(0,0) = 2e\I(Q = 0)|. 

The density of states for the one/double dots sytem is 
even (Eq. (17)) so, according to Eq. (24), the harmon- 
ics of the current have symmetry properties : I(u>) and 
I(3u) are even and I{2u) is odd. 

The harmonics of the current plotted in Figs 7 and 8 
are normalized by 2eirpT\Tj,\ 2 ■ 



A. Single dot problem 

The behavior of the finite frequency higher harmonics 
of the current I(mui) is examined where lu is the driving 
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frequency of the AC-field. 



The stationary current is computed from Eq. (25) 



arctan 



(Ed — Ep) 



arctan 



(E D -ojq- nu>) 



(26) 



Fig. 7 shows the zero-frequency component and the 
first three harmonics as a function of the ratio luv/lu for 
different values of the substrate- molecule coupling which 
defines the effective width of the single energy level of 
dot, irpsT s . The zeroth harmonic represents a general- 
ized Current- Voltage, I-V, curve. For a dot with a sharp 
level structure, the I-V curve displays a staircase struc- 
ture with steps width corresponding to either Aluq = u> or 
Aw = 2w associated with photo- assisted transport. De- 
pending on the overlap between the Fermi distribution 
of the substrate and the band-width of the impurity the 
zeroth frequency harmonic can be either positive or neg- 
ative leading to current inversion at large negative DC 
bias. The Ohmic regime is recovered for large Ts which 
corresponds to the physical limit of a metal (See inset on 
Fig. 7a). 



A general feature of the higher harmonics is the pres- 
ence of maxima and minima whose number increases with 
the order of the harmonic. These oscillations are damped 
as Ts increases and the higher order responses vanish 
while the zeroth harmonic survives. This has important 
implications for the experiment suggested in reference 
[20], because large Tg mimics in our model a metallized 
dot, precisely the situation where Kochanski suggested 
there would be no higher-order signal. This also agrees 
with the results on insulating surfaces and stresses the 
importance of exploring the use of higher harmonics for 
more general moZecwZe-electrode interfaces. 




/ 

LUq/UJ 

FIG. 7. Stationary current and its three lowest harmonics 
as a function of the ratio wo for the single dot system. The 
AC bias is fixed to uj\ = 5.0cj. The Fermi energy and the 
dot level are the same than in previous graph (Ef = —1.4a;, 
E 2 = 0.0). Full line, T s = 0.001; dashed line, T s = 1.0; 
dashed-dotted, T s = 3.0; dotted line, T s = 1000.0. Inset : 
stationary current as a function of loq/lj for small values of 
this ratio (r s = 1000.0). 



B. Two— site molecule 



The double dots current harmonics are obtained from 
Eq. (24) using the appropriate F function (see Eqs. (17) 
and (18)) 

Turning now to the discussion of the current harmonics 
for the 2 dot system (Fig. 8), one first notices -as observed 
in the single dot case- that higher order harmonics oscil- 
late more. For comparison with the previous case, the 
2 dot levels have been chosen symmetrically around the 
position of the single dot level in Fig. 7. 

Starting with sharp levels (poor coupling to the sub- 
strate) , the zeroth harmonics has a monotonous staircase 
structure, ranging from negative to positive value. The 
first harmonic is positive, as in the one dot case, and has a 
broader width due to the spacing between the two levels. 
For the next two harmonics, the oscillatory behaviour is 
more prononced than in the one dot case. As in the sin- 
gle dot case, the harmonics have symmetry properties : 
I(un) is even, I(2ui) is odd, I(3u>) is even. 

Upon increasing the levels widths, the amplitude of the 
oscillations first decreases (up to Ts ~ t). For very large 
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coupling (dotted curve on the figure), we recover exactly 
the one dot case : as explained before, for large Ts, the 
density of states of the double site molecule (Eq. (17)) 
corresponds to the density of states in the single dot case 
for a very small broadening t 2 / ' (71775 r|). 




FIG. 8. Stationary current and its three lowest harmonics 
as a function of the ratio loo/lu for the two site molecule. The 
alternative voltage is fixed to u>i = 5.0u>. The Fermi energy 
E F = -1.4a; and t = 3.0. Full line, T s = 0.001; dashed line, 
T s = 1.0; dashed-dotted, F s = 3.0; dotted line, T s = 1000.0. 



So far, the work on photo-assisted shot noise has fo- 
cused mostly on mesoscopics conductors where the effec- 
tive transmission coefficient has little energy level struc- 
ture [32] . One of the purpose of this work was precisely 
to go beyond this limitation. In fact, we have seen that 
the presence of the molecule levels modifies drastically 
the staircase which is observed in the noise derivative as 
a function of bias. As the coupling to the substrate is in- 
creased, the sharp features which we observed get more 
pronounced. 

This corresponds to a poorly conducting system (in- 
sulator). It is relevant to follow Kochansky's intuition 
when analyzing the behaviour of the current harmonics 
as a function of this coupling. Indeed, the harmonics stay 
sharp in the case of an insulator. To our knowledge, this 
is the first time that this is addressed in a quantitative 
manner, albeit with the help of a simple model. The fact 
that the /(3a;) signal has been measured for insulating 
surfaces makes it reasonable to assume that the experi- 
ment with electrode- moZecwte junctions will soon be done. 
It is on this situation that the predictions of this model 
may be more relevant. 

Future directions of research could include a time- 
dependent, self-consistent treatment of Coulomb interac- 
tions in the molecule. Electronic correlations has recently 
been addressed in References [17,33,34]. Another possi- 
ble extension would be to describe the molecule using ab 
initio method, as used in [35-37] where the effect of an 
AC bias in molecular transport is also addressed. 
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VI. CONCLUSIONS 



In this paper, we have studied the photo-assisted shot 
noise for a simple model of metal- molecule substrate junc- 
tion which could prove to be relevant in alternating cur- 
rent STM experiment. The Keldysh method allowed to 
treat the time-dependent response of this system and in 
particular to express the noise in terms of the dressed 
Green's functions of the molecule plus substrate com- 
pound. Contrary to scattering formulations of trans- 
port, the Keldysh approach allows to take into account 
the molecule occupancy correctly. While the present 
deals with non-interacting electrons throughout, the gen- 
eral approach allows to tackle more complicated systems 
where the molecule would include electron-electron in- 
teractions as they occur in carbon nanotubes or more 
complex nano-objects. 



APPENDIX A: UNPERTURBED GREEN'S 
FUNCTIONS 



1. For two electrodes 



The non-equilibrium Green's functions for a site M 
(M — T, S) into a normal metal electrode are: 



9mm{&) = 2i7rp M /(^) 
9mm(^) = -2mtpm[1 - /(ft)] 

9mm(^) = ™p M = gMM (ft) 



(Al) 



where /(ft) is the Fermi distribution and pu is the 
local density of states on site M near the Fermi level. 
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2. For one dot 



Using the definitions of Keldysh Green's functions 
in terms of creation/annihilation operators of electrons 
in the dot, we find the temporal (homogeneous) non 
equilibrium Green's functions: 

9f D {r)=in Q D e- iE ^ 
9nUr) = -i(l-n° D )e- iE ^ 
9dd(t) = iQ(-T)e~ iE ^ 
g r DD (r) = -iQ{T)e~ iE ° T 



(A2) 

where n D is the number electron in the dot at time 0, 
Ed is the energy and 0(i) is the Heaviside function . 

Using a Fourier transform, the expressions of the 
Green's functions in the frequency space look: 



V 



9 + ooW = 2in%- {n _ ED)2 
g D ~ D (<n) = -2i(l-n D ) 

9 a DD m 

9dd(V) 



(n - E D ) 2 + v 2 



Q — Ed — it] 
1 

Q - E D + it] 



(A3) 



where r\ is an infinitesimal positive constant. 



3. For a linear system of D dots (D > 1) 

A linear chain of D dots is considered. The sites are 
labelled by j = 1..D. The real hopping between two near- 
est neighbor sites is t. In such a system, eigen-energies 
and eigenfunctions of the site j are given by: 



E j =-2tcos(^f T 



Kj sin 



in (w) 



where j, I = 1..D. Eigenfunctions are normalized with 
Kj given by: 



\Kj\ 2 



D 
1=1 



D + l 



(A4) 



The non-equilibrium Green's functions can be derived 
using the previous expressions. In position and energy 
space, that yields: 



g a {hM,E) = Y,\Kn 



n=l 



S + 2tcos(-^ T ) - in 



9 + -(h,l 



^) = ^El^| 2 sin(^L) S in(^ 



x x(E n )S(E - E n ) 



(A5) 



where X (E) is the occupation number of energy level E 
(0 or 1) and r\ is an infinitesimal positive constant. The 
retarded Green's function is the complex conjugate of the 
advanced one. The — h Green's function is obtained by 
replacing X (E n ) by —(1 — \{E n )) in the H — expression. 

In this proposal, only the case D = 2 is studied. The 
expressions of the normalization constants are \Kx\ 2 = 
\K 2 \ 2 — 2/3. The eigenfunctions and eigenvalues for each 
dot are given by: 

( Ei — — t E2 = t 

(A6) 



*i(0 = > /fsin(^) * 2 (/) = V / |sin(M) 
Introducing these expressions in (A5) yields: 



= 3 



sin(2iL)sin(^) sin sin (32a) 



E + t-ir] 



E-t-ir/ 



+- n?\ 4 * 7r 
9h.iS E ) = — 



fWf )M-,WE + n 



+ sin 



2ttZi 



sin 



^ x (t)6(E-t)] (A7) 



For the two dots system, only the Green's functions for 
h, h = 1, 2 are needed: 



1 



1 



Z77 



S?i(£)=.9 2 a 2 (£) = - I 

gUE) = gUE) = - 

yi2\ j U2i\ ) 2 [E + t-ir] E-t-irj] 

g+-(E) = g+-(E) = zir [ X (-t)S(E + t) + X (t)6(E - t)} 
g+-(E) = g+-(E) = ztt [ X (-t)S(E + t) - X (t)S(E - t)} 

(A8) 

APPENDIX B: DRESSED GREEN'S FUNCTIONS 

In this section, dressed Green's functions that enter in 
the expression of the current and noise are derived. In 
the general case of a linear chain of N dots, it is pos- 
sible to compute the dressed Green's functions using a 
linear set of Dyson equations. As explained in the previ- 
ous section, the usual scheme is to isolate the sub-system 
composed by the substrate and the chain and then to ex- 
press the dressed Green's functions of the point contacts 
in term of the unperturbed Green's functions of the iso- 
lated chain and of the substrate. Standard unperturbed 



10 



Green's functions of a metallic electrode are used for the 
substrate (Appendix A 1) including the Fermi energy Ep 
and the density of states ps of the substrate. 



2. For a linear system of D dots (D > 1) 



1. For one dot 

Here, only the system composed by the substrate (S) 
and one dot is considered. The hopping between the sub- 
strate and the molecule will be denoted by T$ and will 
be chosen to be real. The dressed Green's functions have 
homogeneity in time due to the time independence of the 
hopping Hamiltonian. It also permits to use directly the 
Dyson equation in frequency. 

For the advanced and retarded Green's function, using 
two Dyson equations rises: 



l T iTrpsT 2 s9 a o r D m 



(Bl) 



DD\ 



where = Ts the self energy associated to the 

hopping between the substrate and the dot and the 
upper/downer sign stands for the advanced/retarded 
Green's functions. Little g denote the unperturbed non- 
equilibrium Green's functions for an isolated dot. Their 
expressions are given in Appendix A 2. Introducing these 
expressions in the dressed Green's function of the point 
contact gives: 

> (f7) = n-E DT i(r, + *PsT%) m 

where r) is an infinitesimal positive constant. This ex- 
pression show the effect of the substrate. Due to the 
hopping Ts between the dot and the electrode, a broad- 
ening of the energy level of the dot appears. 

The limit when 77 tends to is straightforward because 
these Green's functions have no real poles: 

The same technique allows to obtain the spectral 
dressed Green's functions non-perturbatively: 

9dd (n) + 2t7T Ps rl\ g - DD (n)\ 2 f(n-E F ) 



G+- D (Q) = 



G D + D m 



\l-mpsT%g% D {^ 

-2^ Ps r|| 3 ^(f>)| 2 (i 



f(Q - E F )) 



\l-inpsr%jfo D (tt)\* 



(B4) 



Taking the limit when 77 — > gives: 

/(H - E F ) 



G+- D (n) = 2inpsT 2 i 



S {n _ Ed) 2 + {7rpsT 2 s) 2 

G D +(n) = -2inp s T 2 s 



(1 - m - E F )) 

2 \2 



(n - E D f + (irp s T%) 



(B5) 



In order to compute the dressed Green's functions in 
a non-perturbative way between the substrate and the 
linear chain of dots, the system is studied in the absence 
of the tip. The D dots are labelled by 1..D, where the 
site 1 is the point contact between the substrate and the 
chain and the site D is the point contact between the 
chain and the tip. Green's functions for a linear chain of 
dots are derived in the appendix A 3 and will be denoted 
by a little g. Here, the relevant non-equilibrium Green's 
functions are those for extreme dots. 

In this system, the hopping Hamiltonian between the 
substrate and the chain is time-independent (and so, the 
Green's functions are homogeneous in time), it permits 
to use directly the frequency expression of Dyson equa- 
tion. In the noise and current formulae, only the dressed 
Green's function for the site D is needed. To express in 
a non-pertubative scheme these Green's functions, one 
needs to write a linear set of Dyson equations. Here 
is the expression of the advanced and retarded Keldysh 
Green's functions for the site D: 



= + 9 r 6 a AmV S GT D ^) (B6) 

where = Ts is the self-energy for a single hopping 
and (fl) is a dressed Green's function that has to be 
expand by a new Dyson equation. g T j^(Cl) contains all 
the informations of the isolated linear chain. We find a 
non-perturbative expression for this Green's function to 
all orders of Ts- 



GS£(fi) = 



l±mp s r| 5 [f(0) 



(B7) 



So a formal expression for the non-perturbative 
dressed advanced and retarded Green's functions are ob- 
tained in terms of the unperturbed Green's functions of 
the linear chain of dots: 



Gr.a 
DD 



(fi) 



9iD( n )9D a i(ty 



9DD \ ) T l7r Pst S -, T ■ - F 2 r,a, n , 

l±inp s l s9n ("J 



(B8) 



For the spectral Green's functions, the same technique 
is applied and permits also to have an expression for the 
dressed Green's functions: 
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G+ D (Q)=g+ D (n) 



inpsTl&M l + iirpsTlgl^Q) 
gi a D (») g ^ 1 (»)g 1 + 1 -(») 

+2t7rp s f(n - E F )rl, g r D1 (n)g a 1D (n) 

inp s T s 



l + mp s T 2 s - 



inpsT 2 s9u(ty 

g a nm 



'l-mp s T%gt 1 {n)_ 
G D +(n) = g D + D m 

+tirp s l s 



1 



inp s T 2 s g^{n) l + inpsT 2 s g r u (fl) 

g? D (n)g r D1 (n)9ii + (n) 
: (1 + inpsTlgl^n))^ - inpsTlfoitt)) 

-2inp s (l - f(Sl - E F ))T 2 s g r D1 (n)g a 1D (n) 
tfi(fi) 



'-^^l + inpsTlg^n) 



1 + iirpsTi; 



1 - i7rp s T 2 g^(n) 



(B9) 



In this paper only the case of a double dot molecule 
is considered {D = 2) . The unperturbed Green's func- 
tions for a double dot system are derived in appendix 
A3. Introducing these expressions in Eqs.(B6) and (B9) 
gives: 



G a DD (n) 



fl — iirpsT'g 
~T 2 - Q(n - iTrpsT 2 ] 



G+~ D (n) = 2iwp s T 2 



t 2 f(n - e f ) 
\t 2 -n(n + iTT Ps r 2 )\ 2 



(BIO) 



where t is the hopping between the dots. The dressed 
retarded Green's function is the complex conjugate of 
the advanced one and the — I- Green's function is ob- 
tained by replacing f (fi — Ep) in the -) — expression by 
-(l-f(Q-E F )). 
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